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£S) , Abstract 

{N| ' Within the framework of the standard three neutrino scenario, we derive an exact and 

simple formula of the oscillation probability P(v e — > v^) in constant matter by using a new 
, method. From this formula, it is found that the matter effects can be separated from the 

pure CP violation effects. Furthermore, the oscillation probability can be written in the 
form, P(y e — ► v^) = Acos<5 + _Bsin<5 + C, in the standard parametrization without any 
C^l ' approximation. We also demonstrate that the approximate formula in high-energy can be 



X 



easily reproduced from this as an example. 



o 
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Qh! 1 Introduction 



Just like the quark system, it has been shown from the atmospheric neutrino experiments 
]|] and the solar neutrino experiments Q that neutrinos have finite mass and finite mixing. 
In this situation, it is extremely interesting to investigate the CP phase in the lepton sector. 
Fortunately, recent report from SNO experiment || favors the LMA MSW solutions to the solar 
neutrino problem. This means that the measurements of CP phase may be possible because of 
the large 1-2 mixing angle and the large 1-2 mass difference. 

In order to measure the CP phase, the long-baseline experiments such as the JHF experiment 
Q and the neutrino factory experiments 0] are planned. In the past, the asymmetries APcp = 
P(u a — > up) — P[p a — ► up) and APt = P(u a — > up) — P(up — ► u a ) have been considered as the 
main approach to measure the CP phase 5 |], 0, [8|, || . These are methods to measure the direct CP 
violation term which depends on sin 5. However, the measurement of APqp is not directly related 
to the discovery of CP phase, because of fake CP violating effects from the earth matter. APt is a 
pure T violating observable, but it has its own experimental difficulties. So, alternative approach 
has been recently considered in |ll|, [l^, [l3|, 14]. This is an attempt to obtain the information 



on the CP phase totally from the probabilities itself, not only the direct CP violation term but also 
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the indirect CP violation term which depends on cos 5. In these papers the oscillation probability 
is written approximately in the form, P[y e — ► vS) ~ A cos 5 + B sin 5 + C . The extra information 
which is proportional to cos 5 will lead us to the value of 5 in spite of the matter effect. In 
order to obtain more precise information, it is highly desirable to have an exact expression for 
P{y e —> Vu)- Some attempts to derive the exact formula have been made in the context of three 
neutrino scenarios ||, U, I3> II- These formulae are useful for numerical calculation. However, 



the precise CP dependence of P{y e — ► v^) has not been investigated sufficiently 16]. 

To describe our approach, let us review the work of Naumov |n| and Harrison-Scott ||. The 
Hamiltonian H in matter is related to H in vacuum as 

H = H +^Ldiag(a,0,0), (1) 

where a = 2\^2GpN e E, Gf is Fermi constant and N e is the electron density in matter. In 
particular, taking the products of non-diagonal elements, 

lm(H e ^H, lT H Te ) = lm(H e ^H, lT H Te ), (2) 

one obtains the following identity, which we call Naumov-Harrison-Scott identity, 

A12A23A31 J = A12A23A31 J, (3) 

in CP-odd part, where Ay = mf — mj, J = ImJ^ is Jarlskog factor ppfl , J^g = UaiU^U^Upj 



and U is the Maki-Nakagawa-Sakata (MNS) matrix [21|. Here the quantities expressed by the 



tilde include the matter effects. From this identity, J can be expressed in terms of effective masses 



and the parameters of the Hamiltonian in vacuum. As effective masses shown in ||15| , |16| , 17 1 do 
not depend on the CP phase, J can be completely expressed by the linear term in sin 5. The 
reason why the CP dependence becomes simple is that J is the products of four U's. Complicated 
matter effects included in a U are partially canceled in J. [] 

In this letter, we calculate UU* . We use some matter invariant identities such as Naumov- 
Harrison-Scott identity and express not only J but also Re Jf? with the effective masses and the 
parameters in vacuum. The exact formula obtained in this method is very simple and the matter 
effects come in only through effective masses. We show that Re J£* has only a linear term in cos 5. 
That is, we prove that the oscillation probability in matter can be written in the form 

P(v e -» Vy) = Aco^ b + Bsin5 + C (4) 

without any approximation. 

Another merit of our result is that the exact formula immediately reduces to the well-known 



approximate formulae both in low-energy |7]] and in high-energy |10| , 11 1. We demonstrate that 
the approximate formula in high-energy can be easily reproduced from our exact formula as an 
example. Finally, we numerically calculate the coefficients A, B and C . 



1 The calculation of a U is performed by diagonalizing H in Ref. [ jl7J . 
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2 Exact Formula of the Oscillation Probability 

The flavor and mass eigenstates are related by the MNS matrix U a i in matter, where a = e, fi, r 
is the flavor index, i = 1, 2, 3 is the mass index. The amplitude for u e to transition is given by 

3 A 

A{y e ^v p ) = Y J K^ L V^ (5) 

and the oscillation probability is also given by 

P(ye u„) = \A(v e -> z^)| 2 , (6) 

from the amplitude, where L stands for the baseline length. 

We note from (||) that the amplitude depends only on the products C/ e j{/*j. One of the 
important points in this letter is that these products can be easily calculated from identities 
which we derive below. 

From the unitarity relation and the other two relations, 

H efM = H eiM =p/(2E), (7) 
H eT H Tfl — H efl H TT = H eT H T/1 — H eil H TT = q/(2E) 2 , (8) 

three identities on the products XJ^IJ*^ can be obtained as follows: 



YJjjr^ = J2 u ^u; t = o, (9) 

i=\ i=l 

3 3 

= XXtfeit^p, (10) 
i=l i=l 

cyclic cyclic 

Y, ^kUeiUU = E rfmlU^U;, = q, (11) 
(ijk) (ijk) 

where p and q are constants determined by the parameters in vacuum and the sum is over 
(ijk) = (123), (231), (312). We use the relation U Ti U ej - U ei U Tj = U* k (detU) etc, obtained from 
the formula U' = U^ 1 = ^/(detf/) -1 , where U represents the cofactor matrix. 
Solving the simultaneous equations for the products U e iU*^ we obtain 

ujj- ti = (12) 

where {ijk) takes (123), (231), (312). From the definition jjjj = UeiUfc(UejUtf)* , the exact for- 
mula of the oscillation probability is given by 

cyclic /A..r\ cyclic /A-t\ 

P(u e - „„) = -4 £ ReJg, sin 2 ^ M^Jsin ^ , (13) 
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where the sum is over (ij) = (12), (23), (31) and 



jkjy N 2 A,A J + |g| 2 + Re(^)(A t + A J ) ( 
eM A 4J A 12 A 23 A 3 i 



J = ■ (15) 

A12A23A31 

We find that the matter effects are confined in the effective masses only. We can obtain the 
probability for antineutrinos, v e — > by exchanging a — > —a and 5 — * —5 in A^- and JJ^ of Eq. 

& 

Let us comment on the relation between our result and that of other authors. The second 



identity ( |T0| ) is also given in Ref. 22]. The third identity ( pD is new and play an important role 
in deriving our exact formula. The similar expression to (|i~2] ) is given in Ref. |ll| as the result 
of the calculation of e~ tHL , although the CP phase has not been considered. Next, lm(pq*) in 
(|j~5|) are rewritten as 

Im(p<f) = l/(2Eflm(H ef ,H, lT H re ) = A12A23A31 J, (16) 

from (|7j) and (B). Naumov-Harrison-Scott identity is reproduced by substituting ( |i~6| ) into (15). 

3 Separation of CP odd/even Parts 

In this section, we give a concrete expression for the oscillation probability and then, we study 
the dependence of the oscillation probability on the CP phase. 

First let us consider the constants p and q. We use the standard parametrization 



( C13C12 C13S12 si2,e~ t5 \ 

-C23S12 - S23Sl3Cl2e i5 C23C12 - S23Sl3Sl2e l<5 S23C13 
\ S23S12 - C 2 3Sl3Cl 2 e tS -S23C12 - C23Sl3Si2e tS C23C13 J 



(17) 



where sin^ = Sij, cos 6ij = Cij. In addition, as the neutrino oscillation probabilities do not 

depend on the mass itself, but the mass square differences, we take mf = 0, m\ = A21 and 
ml = A31 without loss of generality. So, p and q are given by 

p = p ie - l5 +p2, q = qie~ i5 + q2, (18) 
where pi and qi are real numbers; 

Pl = (A31 - A 2 l^ 2 )s 2 3Sl3Cl3, P2 = A21S12C12C23C13, (19) 

Ql = -A 3 iA 2 lCi 2 S23Sl3Cl3, 92 = -A31A21S12C12C23C13. (20) 



2 We notice that the expression of UU* in ( Jl^ ) is also derived from Eq. (4) in Ref. after some calculations 
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Then, we have 



\p\ = Pi+P 2 + 2piP2 cos 6, (21) 

kl 2 = gi + g| + 2gig 2 cos5, (22) 

Re(pq*) = piqi + p 2 q 2 + (piQ2 + <HPi) cos 5, (23) 

Im.(pq*) = (p 2 qi -piq 2 )sin5. (24) 

Therefore, the oscillation probability can be written in the form 

P{u e —* Ufj) = A cos 5 + B sin (5 + C, (25) 

from (13)-(|l5|). Note that A, B and C are independent of 8 and the oscillation probability is 

expressed only by the linear terms in cos 5 and sin 5 up to a constant as described below. This is 
one of our main results. Here 



A = E C ^sm 2 (%^), (26) 



cyclic / A • . A 

* = E*-(%r). (27) 

cyclic / A . . r \ 

c = E^ J ( 28 ) 

are expressed by the products of the oscillation part dependent on L and Ay- , B' and Cjj . And 
then, Ay , i?' and Cy are given by 

A _ -4[2pip 2 \j\j + 2q x q 2 + (pig 2 + giggXAj + Aj)] Q 
4J AyAiaAaaAai 

= -2( P2 gi-Pi g2 ) > (3Q) 
A12A23A31 

c ~4[(pj +pl)AjAj + (gj + gf) + (pigi + g 2 p 2 )(Aj + A,)] 

lJ A i7 A 12 A 23 A 3 i 



as the function of the masses and mixing angles. Since the effective masses Aj shown in [15, [0], 17] 
do not depend on 6, the coefficients A, B and C are independent of 5. 

Our analytic result given in (^) should be compared with the result of [|l^] depicted in Fig. 
1 obtained numerically. The trajectory becomes an ellipse in the bi-probability space when 5 
changes from to 2ir. The CP dependence of the exact form of P(u e — ► v^) given in ( |25| ) becomes 
much simpler than the result in jl(]] . By solving (^) for sin 5 and cos 5 one obtain 



sm4 = B(P-C)W+ f -(P ^ (32) 
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Figure 1: An example of CP trajectory We take P for the horizontal axis and P for the 
vertical axis. The value of 5 changes from to 2ir. 

A(P - C) t By/ A* + B*-{P- gg r „v 

COS 5 = ^ 2 + j g2 • ( 33 ) 

Thus, we can determine the value of CP phase except for the ambiguity of the sign from the 
measurement of the probability. The sign ambiguity is understood as follows. If we measure the 
probability of the neutrino at a fixed energy and a baseline, we find the solutions on a "line a" . 
As shown in Fig. 1, there are two intersections X and Y of "line a" with the CP trajectory. This 
is the reason why the ambiguity due to the sign appears in the analytic solutions (]3^ ) and (l33j). 

In order to resolve the sign ambiguity, we need to measure more than two kinds of probabilities, 
for example, neutrino and antineutrino. We denote P and P of the oscillation probabilities for 
neutrino and antineutrino respectively as 

P = AcosS + Bsm5 + C, (34) 
P = Acos5 + Bsm5 + C. (35) 

Then, CP phase can be determined by 



sin 5 



(AP - AP) - {AC - AC) 
AB — AB 



(36) 



cos 5 



{BP - BP) - {BC - BC) 



(37) 



BA-BA 

without the ambiguity of the sign. This means that the solution is at X, the intersect of "line a" 
and "line b" . 

Although the value of CP phase is determined, in principle, in ( |36[ ) and (|37|), there remain 
other ambiguities included in A, B, C and A, B, C. The methods to resolve these ambiguities 



are discussed in the references for example [12, 25, 26, 27 1. We discuss the ambiguities due to the 
sign of mass squared differences in Sec. 5. 
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4 Simple Derivation of Approximate Formula 



In the previous section, we have shown that the exact formula of the oscillation probability 
can be expressed as P{y e — ► vS) = A cos 5 + B sin 5 + C. In this section, we demonstrate that the 
approximate formula seen in |llj is easily derived as an example in the case of mi < m 2 <C m 3 . 
One obtains the approximate formula for other patterns of mass hierarchy in the same way. 

Let us first consider the coefficient B of sin 5. B is expressed in the form of the sum as 
Note that, under the condition x + y + z = 0, the identity 

sin 2x + sin 2y + sin 2z = —4 sin x sin y sin z, (38) 

holds, and B from ( |27| ) is rewritten in the form of product as 

cyclic / a . . r \ 

* - (39) 
(y) v / 

- 

Next, let us consider the coefficient A of cos 5. Under the same condition as in deriving B, the 
identity 

sin 2 x = — (sin x sin y cos z + sin x cos y sin z) . (41) 
holds and A is rewritten as 



cyclic / A • • T \ 

* = E^ sin 2 (42) 



""""Mil 



cyclic 

E 



U + AU) cos ^ sin ^ sin M . (43) 



Substituting @J) and © for p and ? in (H)-©, A, B and C are rewritten with the masses 
and the mixings as 

A = £ -^A 21 [A3 1 A fc (A v A 31 ) + 4 1) ] CQS (^L\ gm f^A gm f^A 

B = J r sin ( sin ( sin ( *g£) , (45) 

A 12 A 23 A 3 i r \ 4E ) \ 4E ) \ 4E )> { ' 

^ A ii A 12 A 23 A 3 i V 4S / ' 
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where J r = si 2 ci 2 s 23 c 23 s 13 c\ 3 , and 

A® = A 21 {A 31 X k (cl 2 -s 2 12 ) + Xts 2 12 -Al 1 c 2 12 }, (47) 

C\f = A 21 A 31 {-\ l {\ j sj 2 + A 31 cj 2 ) - \ 3 (\ lS l 2 + A 31 cl 2 )}sl 3 c 2 13 , (48) 

Cf a) = A| 1 (A iS f 2 + A 3lC ? 2 )(A iS ? 2 + A3ic? 2 ) S2 2 3C?3, (49) 

4 2b) = A2 1 (A,-A 31 )(A,-A3 1 )^ 2 44c?3- (50) 



Note that these expressions are still exact. In the limit of small A 2 i, terms given in Eqs. p7|)-(pOD 
are higher order in A 2 i and can be ignored. The superscripts of A and C stand for the power of 
A21, and (2a) represents the term proportional to s 2 3 and (26) is the term independent of s 2 3 . 

Finally, we obtain the well known approximate formula by neglecting the smallest effective 
mass. In the high energy neutrino the smallest effective mass is Ai ~ A 2l . Other effective masses 
A2 and A3, correspond to a or A31. Accordingly, A, B and C are approximated by 

8J r A 2 iA 3 i . (A 31 L\ . (aL\ . / (A31 -a)L \ 

B = — sin sin sm , (52) 

a(A 31 -a) \ AE J \4E J \ 4E J ' K ' 

4A li , : 2 2 2 2 / (A 3 i- a)L 



C = , A 3 \, 2 slss 2 3 c 2 3 sin 2 ( ^\^ ), (53) 



under the condition A21/A31 < S13. When S13 is smaller than (A21/A31), the term inde- 
pendent of S13 becomes the dominant term. Although the approximate formula derived here is 



in agreement with the ones seen in [10, 11], the derivation is rather simple. Moreover, one can 



easily reproduce the approximate formula in low-energy P]. 

5 Numerical Analysis of CP odd/even Part 

In this section, we investigate the values of the coefficients A, B and C in cases of neutrino 
and antineutrino using the exact expressions. We also investigate them changing the signs of A31 
and A 2 i- 

In this numerical analysis, we take Q\ 2 = vr/4, |A 2 i| = 10 -4 eV 2 , #23 = 7r /4 and IA31I = 
3 x 10~ 3 eV 2 to be consistent with the LMA MSW solution to the solar neutrino problem pi |j 
and the zenith-angle dependences of atmospheric neutrinos fl|. We also take #13 = 0.05 within the 



upper bound of CHOOZ experiment [29|. The baseline length is typically taken to be L = 2900km 
and the matter density is taken to be 3.2g/cm 3 . 

In Fig. 2 we show the coefficients A, B and C changing with the energy E. We observe that 
the sign of A is opposite for example in Fig. 2(a) and (d). We also observe that A and B have 
the opposite sign but C has the same sign comparing Fig. 2(a) with (e). In addition, some peaks 
have appeared in all graphs of Fig. 2 with the change of energy. In case of A31 > 0, the peaks 



S 



10" 



(a) v e -> ^,mi < m 2 < m 3 



4 
3 
2 
1 

-1 
-2 




10" 



2 5 10 20 50 
E (GeV) 

(c) v e -> z/„, m 3 <C mi < m 2 




5 10 20 50 
E (GeV) 



_ 2 (e) v e -> v„, m 2 < mi « m 3 




10" 



2 5 10 20 50 
E (GeV) 

(g) v e -> iy, m 3 <C m 2 < mi 




2 5 10 20 50 
E (GeV) 



10" 



(b) v e -> ^,mi < m 2 < m 3 
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Figure 2: A,B,C at L = 2900 km The graphs of the left and right side correspond to the 
neutrino and the antineutrino respectively. The solid lines, the dotted lines and the dashed lines 
are for A, B and C in all graphs. And from top to bottom, (A 3 i > 0, A 2 i > 0), (A 3 i < 0, A 2 i > 0), 
(A 3 i > 0, A 2i < 0) and (A 3i < 0, A 2 i < 0) cases. 
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around 6GeV in Fig. 2(a) for neutrinos are enhanced compared with those in Fig. 2(b) for 
antineutrinos. Inversely, in case of A31 < 0, the peaks in Fig. 2(d) for antineutrinos are enhanced 
compared with those in Fig. 2(c) for neutrinos. 

These features are understood qualitatively from the approximate formula (|5l|)-(|53|). First let 
us consider the sign of A, B and C . As we found from (|5l|)-(|53|), when the signs of both A31 and 
a change, the sign of A becomes opposite and the signs of B and C do not change. On the other 
hand, when the sign of A21 changes, the signs of both A and B change while the sign of C does 
not change. Next, let us consider the magnitude of the peak around 6 GeV. These are strongly 
affected by the denominator (A31 — a). Since the signs of A31 and a are opposite in Fig. 2(a) 
and (d), the denominator (A31 — a) becomes small and the magnitude of the peaks are enhanced. 
On the other hand, since the signs of A31 and a are the same in Fig. 2(b) and (c), the peaks are 
suppressed. Finally, let us explain the position of the peak in Fig. 2(a) and (d) around 6 GeV. 
Roughly, the peak energy is determined by the following: 



sin 



1.27 



A. 



31 



leV 2 



L 



1 km 



— V 

lGeVy 



1 -> E ~ 6 GeV (at L = 2900 km) . (54) 



As pointed out by Parke and Weiler H, and Lipari [13], the peak energy is lower than the energy 
of 1-3 MSW resonance since the baseline length is short compared with the earth diameter. The 
above discussions on Fig. 2(a)-(e) can be applied to other figures. 

We have studied how the magnitude of A, B and C change due to the sign of the mass squared 



differences. In the case of mi < m<i <C 7713, the coefficients have been investigated in Ref. [1C] 
by using the approximate formula. These correspond to Fig. 2(a) and (b). The sign of A31 is 
determined from the leading term C as pointed out by many authors (for example see [pp . On 
the other hand, the sign of A21 is determined from next leading terms A or B. This means that 
the sign of A21 is simultaneously determined in addition to the CP phase. It may be interesting 
as the first observation of the sign of A21 using artificial neutrino beam. 



6 Summary 

We have studied neutrino oscillations in constant matter within the framework of the three 
neutrino scenario. We summarize the results obtained in this letter. 

(i) We have derived an exact expression of the oscillation probability by using a new method. 
We have calculated UU* from the identities without directly calculating single U. Not only 
the derivation but also the result becomes simple and the matter effects enter only through 
effective masses. 

(ii) We have obtained the CP dependence of the oscillation probability exactly by using the 
standard parametrization. It has been shown that the oscillation probability is in the form, 



Piye 



A cos 8 + B sin 5 + C. We have also demonstrated that the approximate 



formula in high-energy can be easily reproduced from our result. 
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Finally, let us comment on the oscillation probabilities for other channels. These probabilities 
are easily derived in the same way as P(u e — ► v^). P{y e — > v T ) has the same CP dependence 
as P{y e — * Vu). However, P(vn — * v T ) has the term which depends on cos2<5 in addition to the 
linear terms in sin 5 and cos 5. 
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